Uniformly convergent non-standard finite difference methods for self-adjoint singular perturbation problems
نویسندگان
چکیده
We design non-standard finite difference schemes for self-adjoint singularly perturbed two-point boundary value problems. Essential physical properties (e.g., dissipativity) of the solutions of such problems are captured in the schemes by an appropriate renormalization of the denominator of the discrete derivative. The schemes are analyzed for -uniform convergence. Several numerical examples are given to support the predicted theory. © 2005 Elsevier B.V. All rights reserved. MSC: 65L10; 65L12; 65L70; 65L99
منابع مشابه
Wavelet optimized finite difference method using interpolating wavelets for self-adjoint singularly perturbed problems
We design a wavelet optimized finite difference (WOFD) scheme for solving self-adjoint singularly perturbed boundary value problems. The method is based on an interpolating wavelet transform using polynomial interpolation on dyadic grids. Small dissipation of the solution is captured significantly using an adaptive grid. The adaptive feature is performed automatically by thresholding the wavele...
متن کاملHigh order fitted operator numerical method for self-adjoint singular perturbation problems
We consider self-adjoint singularly perturbed two-point boundary value problems in conservation form. Highest possible order of uniform convergence for such problems achieved hitherto, via fitted operator methods, was one (see, e.g., [Doolan et al. Uniform numerical methods for problems with initial and boundary layers, Boole Press, Dublin, 1980], p. 121]). Reducing the original problem into th...
متن کاملOn Approximate Stationary Radial Solutions for a Class of Boundary Value Problems Arising in Epitaxial Growth Theory
In this paper, we consider a non-self-adjoint, singular, nonlinear fourth order boundary value problem which arises in the theory of epitaxial growth. It is possible to reduce the fourth order equation to a singular boundary value problem of second order given by w''-1/r w'=w^2/(2r^2 )+1/2 λ r^2. The problem depends on the parameter λ and admits multiple solutions. Therefore, it is difficult to...
متن کاملReliable Finite Element Methods for Self-adjoint Singular Perturbation Problems
It is well known that the standard finite element method based on the space Vh of continuous piecewise linear functions is not reliable in solving singular perturbation problems. It is also known that the solution of a two-point boundary-value singular perturbation problem admits a decomposition into a regular part and a finite linear combination of explicit singular functions. Taking into acco...
متن کاملParameter-uniform Fitted Operator B-spline Collocation Method for Self-adjoint Singularly Perturbed Two-point Boundary Value Problems
Under these conditions the operator L admits the maximum principle [1]. These types of problems, in which a small parameter multiplies the highest derivative, are known as singular perturbation problems. They arise in the mathematical modeling of physical and chemical processes, for instance, reaction diffusion processes, chemical reactor theory, fluid mechanics, quantum mechanics, fluid dynami...
متن کامل